Landau-Zener transitions in the presence of 

harmonic noise 

Matthias Kraft 1 ' 2 , Stephan Burkhardt 1 , Riccardo Mannella 3 , and 

Sandro Wimberger 1 

Institut fur theoretische Physik, Universitdt Heidelberg, Philosophenweg 19, 69120 Heidelberg, 

Germany 

2 Blackett Laboratory, Imperial College London, South Kensington Campus, SW7 2AZ London, 

UK 

3 Dipartimento di Fisica E. Fermi, Universita di Pisa, Largo Pontecorvo 3, 56127 Pisa, Italy 

March 13, 2013 



We study the influence of off-diagonal harmonic noise on transitions in a 
Landau-Zener model. We demonstrate that the harmonic noise can change 
the transition probabilities substantially and that its impact depends strongly 
on the characteristic frequency of the noise. In the underdamped regime of 
the noise process, its effect is compared with the one of a deterministic si- 
nusoidally oscillating function. While altering the properties of the noise 
process allows one to engineer the transitions probabilities, driving the sys- 
tem with a deterministic sinusoidal function can result in larger and more 
controlled changes of the transition probability. This may be relevant for 
realistic implementations of our model with Bose-Einstein condensates in 
noise-driven optical lattices. 

1 Introduction 

The Landau-Zener model is one of the simplest, yet non-trivial time-dependent problems 
in quantum mechanics. Originally solved in 1932 [HH], it has since provided a valuable 
tool in the study of two state quantum systems [5 , 6j . In recent years it experienced a 
renaissance amongst both, theorists and experimentalists. On the experimental side, new 
techniques allowed to test theoretical predictions with improved precision [Tltll] . On the 
theoretical side, in-depth studies on finite time effects in Landau-Zener transitions |12pi3j 
and on "noisy" Landau-Zener transitions [T1HT6] were conducted. 
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Applications of the Landau-Zener model are versatile; it can provide a tool to model 
and manipulate "qubits" in circuit QED [17] . to calculate transitions in semiconductor 
quantum dots |181I19| or enhance the experimental control over Bose-Einstein conden- 
sates in optical lattices [20]. Along those lines, substantial effort has been made to find 
new ways to control Landau-Zener transitions by, for example, developing new experi- 
mental protocols [21J. Another approach to control the transitions, which we will also 
follow in this letter, is the introduction of noise into the model [22J. A different approach 
to the role of noise in bistable quantum system can be found in |23j. 

In this letter, we present a modified version of the Landau-Zener model containing 
harmonic noise that directly originates from studying the behavior of Bose-Einstein 
condensates in noise-driven optical lattices (see [23] as well as [25]). We demonstrate 
that in this model, the transition rates can be effectively controlled by changing the noise 
properties and that this control can be even enhanced by replacing the noise process by 
a deterministically oscillating phase. 

This letter is organized as follows: first we present our modified Landau-Zener model, 
then we introduce harmonic noise and discuss its most important properties. In section 
[3 we give numerical results for the diabatic transition probabilities. We conclude with 
further perspectives on experimental realizations in section 21 

2 "Noisy" Landau-Zener model and Harmonic Noise 
2.1 Landau-Zener model with off-diagonal noise 

In this section, we will briefly introduce our modified "noisy" Landau-Zener model. The 
Hamiltonian of our model, motivated by studies of BECs in optical lattices, reads 



Here, the main difference to the standard Landau-Zener model [IH3] is the appearance 
of a complex phase on the off-diagonal. As <p represents our noise variable (more details 
in section [272]) . this is also the term responsible for the stochastic nature of our model. 
The main difference to commonly studied "noisy" Landau-Zener models [14[I15] is that 
the coupling term on the off-diagonal is composed of a constant, as well as a noise- 
driven term. Physically, this model can be used to analyze Bose-Einstein condensates in 
tilted optical lattices. In this case, each of the two terms on the off-diagonal represents a 
sinusoidal optical lattice. While the constant term stands for a static lattice, the complex 
phase term stands for a second optical lattice of the same wavelength [24J , shifted by a 
phase cp. 

The states |0) d = (1,0) and |l) d = (0,1)' are called the diabatic basis of the 
system (eigenstates of JI]) with Vo = 0), while the time-dependent energy eigenstates 
\0}a(t),\l) a (t) are called the adiabatic basis of the system [13]. The diabatic and the 
adiabatic eigenstates are identical in the limit t = ±oo. However, around t = the 
adiabatic eigenstates go through an avoided crossing, whereas the diabatic eigenstates 




(1) 
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cross each other. Thus the correspondence between the two bases "switches" such that 
|l) a (t = — oo) = whereas \l) a (t = +00) = \0)d- The instantaneous stochastic eigen- 
values (adiabatic states) of this model can be found by diagonalizing the Hamiltonian 
flU ; we then define an effective band gap0 at t = by averaging over the stochastic phase 
[M]. This yields, 

A£ e ff = W2(cos(<K0)) + l)>, (2) 

where (...) denotes an average over the equilibrium distribution of 4> (see fll])). 

The most interesting property of this system are the diabatic and adiabatic transition 
probabilities. The diabatic transition probability (Pd.tra) is the probability to end up in 
\l)d at t = +00 if the system was initially prepared in the opposite state \0}d at t = —00, 
or vice versa. The adiabatic transition probabilities accordingly refer to changes of the 
adiabatic state. Since the diabatic and adiabatic eigenstates are identical at t = —00 but 
are "swapped" at t = +00, the adiabatic transition probability is just P a ,tra = 1 — -Pd,tra- 
These transition probabilities can be computed exactly for the original Landau-Zener 
problem (<p(t) = in equation (JT])) [1-4J. The diabatic transition probability reads 
-Pd,tra = 1 — exp(— with Vq being the band gap at t = 0. Using the bandgap given 
in ([2]), the formula can also give an estimate of the diabatic transition probability for 

the Hamiltonian ([1]) by replacing Vq with AE^ S , i.e. P est = 1 — exp(— n ^p° s )■ 
2.2 Harmonic Noise 

Previous studies of noise-driven Landau-Zener problems were mainly concerned with the 
influence of white [2] or exponentially correlated noise |15}I16]. To better understand 
the interplay between noise and the timescales of the system itself, in this letter the 
effect of a noise process with a characteristic frequency is investigated. Here, we will 
give a short overview over this so-called harmonic noise process. 

The Langevin equation defining harmonic noise can be written as [261 127] 

d> =v 

r (3) 

v = - 2Tv - ul<$> + V4TT£(t). 

From these equations, we can see that the harmonic noise process is nothing but a 
damped harmonic oscillator driven by a white noise process £(t). The properties of 
this white noise process are determined by (£(£)) = and = S(t — t'). The 

equilibrium distributions of the two noise variables 4> and v are independent and of a 
gaussian shape with [26j 

W)) = (^f) = ^2 W 

("(*)> = (Ht) 2 ) = T. (5) 



1 The band gap is defined as the minimal difference between the two instantaneous eigenvalues of 
Hamiltonian |l} and thus a function of the stochastic phase <j>(t). 
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Figure 1: Power spectrum of a harmonic 
noise process in the under- 
damped (red solid line) and 
over-damped case (dotted blue 
line). Noise parameters are: 
T = 1.0, T = 0.5, w = 5.0 
(under-damped case) and T = 
0.01, T = 2.5, w = 1.0 (over- 
damped case). 

At first it might seem counterintuitive that these distributions do not depend on the 
damping coefficient T. The reason for this is the form of the noise term in ([3]), which is 
scaled with VT and thus compensates the damping. 

Many properties of the harmonic noise process can be understood through the analogy 
to the damped harmonic oscillator. Just as its non-stochastic counterpart, the behav- 
ior of the harmonic noise can be separated into two different regimes: the oscillating 
(under-damped) and non-oscillating (over-damped) one. If Wq > 2T 2 , the process shows 
oscillations, while for the opposite case, its behavior is closer to an exponentially corre- 
lated noise process. 

The difference between these two regimes becomes evident when looking at the power 
spectrum of the noise variable <f>(t). This power spectrum is given by 

2TT 

Tr{4r z Lu z + (or — lUq) r) 

and can be seen in figure [TJ For the oscillating case, we can see that this response 
function does have a clear maximum at oj\ = v/wq — ^r 2 , while for the non-oscillating 
case, the maximum is assumed at cj = 0. For small damping coefficients, the width of 
the maximum is proportional to 2T. This power spectrum is very similar to the response 
function of a damped harmonic oscillator driven with the frequency oj. 

Without damping (r = 0), Q turns into the differential equation of a harmonic 
oscillator and the power spectrum (0) correspondingly turns into a delta distribution 
5(uj — u)q). If the damping time 1/T is very small compared to the timescale of the 
system, the noise process is approximated by a simple harmonic oscillation. In order 
to conserve the total power (integral over the power spectrum) of the process, this 
approximation should take the following form |28| : 

/ j 1 

4>(t)det = 2W — sin(wo* + <fo) = 2 V / Var(0) sin(w t + 4>o)- (7) 
V w o 
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3 Numerical Results 



In this section, we present numerical results for the diabatic transition probability in 
dependence of the characteristic frequency of harmonic noise. The results have been ob- 
tained by numerically integrating the time-dependent Schrodinger equation correspond- 
ing to the Hamiltonian ([1]); the initial state was chosen to be the diabatic ground state 
of the model. The integration time was long compared to the transition time |12lll3j 
such that finite time effects were minimal. We averaged the diabatic transition proba- 
bility over the last ten percent of the integration time to estimate its asymptotic value 
at t = oo. For details see [28J. In the following "transition probability" will always refer 
to the diabatic transition probability. 
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Figure 2: Diabatic transition probability versus the dimensionless noise frequency ojq. 

The top figure shows results for Var(0) = 0.1 (red dotted, triangles), Var(0) = 
0.5 (green dashed, crosses) and Var(^) = 2 (blue solid, lines). The other 
parameters are Fq = 0.5, V = 0.05, and Vq = 0.75. The vertical solid black line 
at coq = 0.07 marks the point after which the noise shows oscillatory behavior. 
The horizontal lines give the saturation value around which the transition 
probability oscillates. 

The bottom figure shows results for T = 0.0 (red dotted, triangles), F = 0.05 
(green dashed, crosses) and T = 2.5 (blue solid, lines). The other parameters 
are Fq = 0.5,Var((/>) = 0.5, and Vq = 0.75. The vertical lines (solid green at 
ojq = 0.07 for T = 0.05 and solid blue at dq = 3.53 for T = 2.5) mark the 
points after which the noise shows oscillatory behavior. The horizontal lines 
again give the saturation value of the transition probability. 
The error bars are only shown for every 8th point and give the error of the 
expectation value. 
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Figure [2] shows the transition probability versus the dimensionless noise frequency 
uq for different variances (top) and damping coefficients V (bottom). The transition 
probability strongly depends on the noise frequency. Comparison with the transition 
probability P cst ~ 0.989 — 0.999 (for all graphs) estimated using the effective band gap 
from ([2]) , reveals that harmonic noise strongly reduces the diabatic transition probability, 
i.e. it increases the adiabatic transition probability (P a ,tra = 1 — -Fd.tra) to the upper 
band. This opens up the possibility to influence the transition probability by changing 
the noise parameters (see also [22]). 

Taking a more detailed look at figure [2j we see that for all six data sets the influence of 
the noise is small for small frequencies, but the transition probabilities strongly decrease 
between uj$ ~ 0.5 — 1.5 leading to clear local minima (apart from the data with V = 2.5). 
The positions of those minima depend linearly on the effective band gap; this effect is 
studied in detail in |24j . 

Disregarding the data for V = 2.5, the transition probability oscillates around a satu- 
ration value for increasing loq. Figure [2] (top) shows that this saturation value depends 
on the variance of the noise process, i.e. on the noise strength. For small variances 
the noise influence is relatively small; increasing the variance enhances the noise influ- 
ence, until the saturation value drops to m 1/2 for Var(^>) = 2. Strong noise thus leads 
to an incoherent superposition of the two energy states, a fact already discovered by 
Kayanuma [14^115] for similar LZ models. This supports the intuitive notion that strong 
noise leads to a complete mixing of the energy levels where the band gap is minimal. 

Figure [2] (bottom) shows that changing the damping coefficient T does not change the 
saturation value as much. However, it can be seen that changing the damping coefficient 
controls the sharpness and height of the oscillations around the saturation value (they 
are completely damped out for T = 2.5). 

Analyzing the two extremes, r = ((f) describes sinusoidal motion and we average 
over amplitude and phase according to Q) and T = 2.5, we can see that the transition 
probability for T = still shows the same characteristic decrease and oscillations. The 
transition probability for T = 2.5 does not show those oscillations and the initial decrease 
happens at a slower rate. Nevertheless, the saturation value is almost the same as for 
T = 0.05. We can thus conclude that the existence of this saturation value is a general 
feature of harmonic noise, independent of the noise parameters we are working with. On 
the other hand, the rapid initial decrease in the transition probability and the oscillations 
around a saturation value are due to the sinusoidal nature of harmonic noise. This leads 
us to the next figure, in which we compare harmonic noise with a deterministically 
oscillating phase. 

Figure [3] shows the transition probability versus the dimensionless frequency ujq for a 
harmonic noise process and a deterministically oscillating phase as given in Eq. ([7]) (we 
give both an average over </>o and a single realization with <j)Q = 0). It can be seen that 
there is a good qualitative agreement between the data obtained with harmonic noise and 
with a deterministic phase. The LZ model with a deterministic phase nicely reproduces 
the initial decrease in the transition probability and the frequency of the oscillations 
around the saturation value. However, the oscillations are much more prominent for 
the deterministic phase and the positions of the first few minima/maxima seem slightly 
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Figure 3: Diabatic transition probability versus the dimensionless noise frequency 
ojq. We show data for the "noisy" LZ-model with harmonic noise (r = 
O.O5,Var(0) = 0.5) (green dashed line, crosses), with a deterministic phase 
4> = v2 sin(woi + (fro) averaged over </>o(red dotted line, triangles) and a single 
realization of the deterministic phase noise with (fro = 0.0 (blue solid line). The 
horizontal lines give the saturation value (dashed green line for the noisy LZ 
model with harmonic noise, red dotted line for the deterministic phase and 
blue solid line for the single realization of the deterministic phase). For all 
three, Vo = 0.75 and F$ = 0.5. (fro has been chosen from a flat distribution 
between and 2ir. We averaged over 100 realizations. Error bars as in figGJ 



shifted to lower values of the frequency. 

A possible explanation for this shift and why it vanishes at high frequencies could be 
that the noise does not behave oscillatory at those minima yet and that our model thus 
breaks down. However, the noise already shows oscillatory behaviour for ojq > y/2T = 
0.07, so it will clearly be highly oscillatory at frequencies where the minima/maxima 
occur (first minimum is at ojq ~ 1.6). One will reach the same conclusion by looking at 
the power spectrum of the harmonic noise © . For the parameters in figure [3] and with 
wo = 1.6, the power spectrum is very narrowly peaked at ui\ = \J — 2F 2 = 1.5992 and 
almost completely symmetric, thus comparing very well with the "delta peaked power" 
spectrum of the deterministic phase. This can, therefore, not be the reason for the 
discrepancies. 

The real shortcoming of our model is a different one, namely that the amplitude of 
harmonic noise varies, whereas the one of the deterministic phase stays constant. Prelim- 
inary numerical calculations have shown that different amplitudes of the deterministic 
phase lead to slightly different transition probabilities, even though the overall oscillatory 
behaviour and the approximate position of the minima/maxima stays intact. Harmonic 
noise should, in the highly oscillatory regime, thus be more accurately described by a 
harmonic oscillator with a randomly distributed amplitude that follows a gaussian distri- 
bution. However, for the sake of simplicity of our model, we refrained from introducing 
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a stochastic amplitude. 

Amplitude and sharpness of the oscillations render the deterministic phase model even 
more suited to control the transition probability in the "noisy" LZ model than harmonic 
noise. Focussing on the transition probability around ojq ~ 3, we see that a slight change 
in the frequency results in a huge amplification (reduction) of the transition probability. 
In numbers, the transition probability rises from 0.153 at ujq = 3.12 to 0.958 at loq = 3.39 
for the deterministic phase model averaged over the random phase <pQ. This effect will be 
greatly enhanced if one is furthermore able to control the random phase (pQ. In this case, 
the transition probability rises from 0.005 at ujq = 3.12 to 0.997 at uq = 3.39 (about 200 
times the value at ojq = 3.12). Even though the difference in the transition probability 
between subsequent minima and maxima is not always as prominent, it is still obvious 
that one can adjust the transition probability over orders of magnitude by only slightly 
changing the frequency of the deterministic phase. It should be noted that we obtained 
these graphs without any optimization to find the "best" parameters that maximize the 
impact of the harmonic noise or the deterministic phase. 

4 Conclusions 

In this letter, we have shown that the here proposed Landau-Zener type model allows 
for greater control of the transition probabilities between the two energy states when 
compared to the standard Landau-Zener model [IHl]- Driving the system with harmonic 
noise, we were able to change the transition probabilities while keeping the coupling 
term Vo as well as the driving Fq constant and only changing the properties of the 
noise. Compared to previous studies of noise-driven Landau-Zener problems [HHTB] . we 
were able to observe a wide range of behavior, as seen in figure [2j We demonstrated 
that the effects most relevant for controlling the transitions, such as the dependence of 
the transition probability on noise frequency, are due to the sinusoidal nature of the 
harmonic noise process (see figure 

While the system under the influence of harmonic noise allows to change the tran- 
sition probability by tuning the properties of the noise process, this can be done even 
more efficiently by using a deterministically oscillating phase (see Eq. ©). Our numer- 
ical simulations revealed that using this deterministic phase, even small changes in the 
frequency cause the diabatic transition probability to change over orders of magnitude. 
This, of course, provides the chance to develop new experimental protocols (see also |21j ) 
that allow for an excellent control over the system parameters. An ideal playground to 
implement such protocols is given by optical realizations [Tllllj or Bose-Einstein conden- 
sates in optical lattices p fTUll2"0H2"T] . 
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